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Abstract 

We discuss the properties of the position space fermion propagator in three dimensional QED 
which has been found previously based on Ward-Takahashi-identity for soft-photon emission ver- 
tex and spectral representation in quenched approximation. There is a new type of mass singularity 
which governs the long distance behaviour.lt leads the propagator vanish at large distance. This 
term corresponds to dynamical mass in position space and expressed as superposition of the propa- 
gator with different mass in momentum space. Our model shows confining property and dynamical 
mass generation for arbitrary coupling constant. Since we used dispersion relation in deriving spec- 
tral function there is a physical mass which sets a mass scale. Low energy behaviour of the 
proagator is modified to decrease by position dependent mass. In the limit of zero infrared cut-off 
the propagator vanishes with new kind of infrared behaviour. 
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I. INTRODUCTION 

Infrared behaviour of the propagator in the presence of massless particle (i.e. photon and 
graviton) had been discussed assuming spectral representation of the propagator and the 
use of Ward-Takahashi-identity[l].It is pointed out that the spectral representation is model 
independent and the asymptotic form of the scalar propagator in QED and Gravity theory 
were determined[l].It has long been thought that the three dimensional QED is a confining 
theory since its infrared divergences is severe [3]. However it has not been known what we 
can observe as a confinement. One of these examples is a vanishment of the propagator in 
the infrared region. In the previous paper we find that the scalar and spinor propagator have 
a new type of mass singularity in the same approach in ref[l] to QED 3 [2]. There are two 
kinds of gauge invariant mass singularities in the evaluation of the 0(e 2 ) spectral function 
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by LSZ reduction formula.After exponentiation of 0(e 2 ) spectral function we get the full 
propagator in position space which is a product of free propagator with physical mass and 
quantum correction. In our model quantum corrections are Coulomb energy and position 
dependent mass which are both logarithmically divergent at long distance. Here confinement 
means that the propagator damps faster than the free propagator with physical mass for 
arbitrary coupling constant. First we study the structure of the propagator with finite in- 
frared cut-off pL.lt is assumed that the three dimensional analysis leads the leading order of 
high temperature expansion results in four dimension [3]. Using Laplace transformation of the 
position dependent mass numerically, we have the integral reresentation of the propagator in 
momentum space as in dispersion theory. As a result we find in momentum space that the 
propagator vanishes in the limit of zero infrared cut-off fi 2 = m 2 —p 2 near p 2 = m 2 .In section 
II Bloch-Nordsieck approximation in three dimension is reviewed. In section III mass singu- 
larity in four and three dimension is compared and numerical analysis in momentum space is 
given. Section IV is for non-perturbative effects as renormalization constant, bare mass which 
are defined in the high-energy limit, and the vacuum expectation value of the composite op- 
erator i^tf) .Renormalization constant Z^ 1 vanishes for arbitrary coupling constant. There 



II. WARD-TAKAHASHI-IDENTITY AND THE SPECTRAL FUNCTION 
A. Position space propagator 

First we consider about the charged particles which emitt and absorb massless pho- 
tons. Usually this process was described by spectral function;transition probability of par- 
ticle into particle and photon state. This method is model independent and helpful in any 
dimension and leads to the so called Bloch-Nordsieck approximation. Multi-photon emitted 
from external line is introduced by ladder type diagrams which satisfy Ward-Takahashi- 
identity.Let us begin by dispersion theoretic description of the propagator[l] 



seems to be a critical coupling constant which separates the phases of 





(1) 
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-QS(p)= I du 2 5(p 2 - c 2 )[ 7 • p Pl (u 2 ) + p 2 {u 2 )} 



IT 



= 1 -PPiip) +P2(p), (2) 

s F (x) = (n\Tip(x)^(o)\n) . (3) 

The field ip is renormalized and is taken to be a spinor with mass m.Here we introduce 
intermediate states that contribute to the spectral function 

p(p) = 7 -PPi(p) +P2(p) 

= (2tt) 2 $ 3 (P ~ Pn) f d 3 x exp(-ip ■ x) (Q\ip(x)\N) (JV|^(0) . (4) 

Total three- momentum of the state \N) is p^.The only intermediates N contain one spinor 
and an arbitrary number of photons. Setting 

\N) = \r;k u ...,k n ), (5) 

where r is the momentum of the spinor r 2 = m 2 ,and hi is the momentum of ith soft photon, we 
have 

p(p 2 ) = J ^ E k\ x ( / w? d ^ )5{k2) '• ••>> 

x (Q|V(x)|r; fci, .., k n ) (r; k u ..k n \^(0)\Q) . (6) 
Here the notations 

n 

(/(*))o=i, (/(*))» =n/(*i) (7) 

i=i 

have been introduced to denote the phase space integral of each photon. Intial sum over e is 
a sum over polarization of photon. To evaluate the contribution of soft-photons,we consider 
when only nth photon is soft. Here we define following matrix element 

T n =(Q\^\r;k u ..k n ). (8) 

We consider T n for k 2 ^ 0,we continue off the photon mass-shell by Lehmann-Symanzik- 
Zimmermann(LSZ)formula: 

T n = e^Tnn, (9) 
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u 



= (n\T^(k)\r;k u ...k n ^J (10) 

lim ( d z xe*v{ik n -x)%(tt\Ti)e»Al' in {x)\r-k l ,...k n _ 1 ) (11) 

J d 3 xexp(ik n -x)(D x + fx 2 ) (n\T^A li (x)\r;k 1 ,...k n - 1 ) (12) 

d 3 xexp(ik n ■ x) {V^T^e^r; k u ../c n _i) , (13) 



% 



provided 

(D x + fi 2 )T^A,(x) = Tip(H x + fi 2 )A,(x) = m-Ux) + (1 - \)d?(d ■ A(x))), 

2 2 

A(D X + !L)( d .A) = 0, + /i 2 )^9(9 • A{x)) = (A - 1)^(0 • (14) 

where the electromagnetic current is 

f(x) = -e^(x)^(x), (15) 

and 1/A is a gauge fixing parameter and /i is a photon mass. In deriving the spectral rep- 
resentation for photon propagator,the Stueckelberg approach was used [1] and the realistic 
condition for transverse component of incoming field is defined 

A T / n (x) = Af(x) + - 2 d,d ■ A m (x), (16) 
/i 

and assumed 

(Tin A 



M^t^-oo - VHK' 111 ~ z ~ d » d ■ Am ^)\- ( 17 ) 

Condition for matrix element in particular is 

(0|^(x)|l> = v^[<0|Af-(x)|l> - z^d, (0\d ■ A m {x)\l)l (18) 

where z is a renormalization constant for londitudinal part [1]. An alternative way to derive 
the reduction formula is an imposing Lorentz condition for the physical state 

d- A ( - +) \phys>=0. (19) 

From the definition (9), (10) T n is seen to satisfy Ward-Takahashi-identity: 

k nfl T%(r, k u ..k n ) = eT n _i(r, h, ../c„_i),r 2 = m 2 , (20) 
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provided by equal-time commutation relations 

^T(iPj,(x)) = -e^(x), 

d;T@j,(x)) = e^(x). (21) 

In the Bolch-Nordsieck approximation we have most singular contribution of photons which 
are emitted from external lines. In perturbation theory one photon matrix element is given 



7i = ( in\T(ip in (x),ie J d x^ in (y)-f^ in (y)A^ n (y))\r; k in 

d 3 yd 3 zS F (x - y)7 M 5 (3) (y - z) exp(i(k -y + r- z))e ll (k ) X)U(r, s) 
(r + k) • 7 + m 



it 



= —le- 



-7 M e M (A;, A) exp(i(A; + r) • x)U(r, s), (22) 



(r + k) 2 — m? 

where U(r,s) is a four-component free particle spinor with positive energy. U(r, s) satisfies 



the relations 



(7 • r — m)U(r, s) = 0, 



^U(r,s)U(r,8) = l^±^. (23) 

In this case the Ward-Takahashi-identity follows 

K T \ = ~ ie — I rr\ (7 • k)U(r, s) 

7 • (r + k) — m 

= -ieU(r,s)=eT , (24) 

provided lowest-order Ward-identity 

7 • k = (7 ■ (r + k) — m) — (7 • k — m). (25) 

For general T n low-energy thorem determines the structure of non-singular term in k n by 
the requirement of gauge invariance of total transion amplitudes under the shift e M — > e M + 
c/c^.They give finite correction to the pole terms for examples in the Bremuthstralung or 
Compton scattering. Detailed discussions are given in ref[l] and non-pole terms are irerevant 
for the single paticle singularity in four- dimension. Under the same assumption in three- 
dimension we have 

T n \ k 2 =0 = e ]—7T~\ Tn-i- (26) 

7 • (r + k n ) - m 



106 



From this relation the n-photon matrix element 



m(x)\r; h, ., k n ) (r; , h,., k n \4>(0)\n) 



(27) 



reduces to the products of lowest-order one photon matrix element 

n 

T n T n = \[T l {k j )T + {k ] ) 1 ,. 



(28) 



In this case the spectral function p(p) in (6) is given by exponentiation of one-photon matrix 
element, which yields a infinite ladder approximation for the propagator. In this way the 
spectral function is given in the followings 

md 2 r 



/Tfld T 
— -p-^exp(ir -x)exp(F), 

F= W(*)M>(r;A#(0)|fi) 

one photon 

/d 3 k _ 
——^5(k 2 )9(k ) exp(ik ■ x) ^ T X T X . 
[Z7T) x,s 

To determine F first we take the trace of T{T\ for simplicity in the infrared. 

d 3 k 



exp(ik ■ x)5(k 2 )6(k ) 



e 

x — tr 
4 



(r + k) • 7 ^ r -7 + m v (r + k) ■ r 
(r + k) 2 — mP 2m ^ (r + k) 2 — m 2 ^ 



Here is the polarization sum 



and the free photon propagator is 



We get 



F = -e 2 



exp(iA;-a;)^(A;o)[5(A; 2 )(- 



m 



5(k 2 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(2tt) 3 ^ v / v u/l v /v (r- k) 2 r-k> v A; 2 J ' 

The second term S(k 2 )/k 2 equals to —5'(k 2 ).Om calculation is the same with the evaluation 
of the imaginay part of the photon propagator. To aviod the infrared divergences which 
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arises in the phase space integral,we must introduce small photon mass p as an infrared 
cut-off. Therefore (28) is modified to 

d 3 k 



/d k 
-^-^ exp(ik ■ x)0(ko) 



x W - ^>((^ + - <« - D^(* 2 - (35) 

Here we assume pi(u 2 )/m = P2{^ 2 ) = p which is valid in the infrared(i.e.7 • r = m).In 
general case there are two kinds of spectral function which is given in the appendix.lt is 
helpful to use function D + (x) to determine F 



D + {x) = [ exp(ik ■ x)d 3 k9(k°)5(k 2 - p 2 ) 

(2ixYi J 



nkdk 

J (k \x\)- 



exp(— p 


x\) 


8ni 


X 





\x\ = V—x 2 . (36) 

If we use parameter trick 

,. f°° , . x / \ \ iexp(ik-x) . . 

hm / daexphik + le) ■ (x + or = ; , (37) 

Jo k-r 

lim / ada exp(i(k + ie) ■ (x + ar)) = . r ,\ (38) 

e^o J - (k-r) 2 

the function F is written in the following form 



/•oo rcc q 

F = ie 2 m 2 / adaD + (x + ar, p) — e 2 / daD + (x + ar, p) — ie 2 (d — l)——D + (x, p) 
Jo Jo dp 2 

_e 2 m 2 / exp(-/i|x|) +| ^ |Ei ^^ |)) _^ 2 _ E .^ |a;|) + (rf _ i) e^^^^^^ (gg) 



87rr 2 p 8nVr^ 8np 

F 1 + F 2 + F g , (40) 



where the function Ei(p \x\) is defined 



Ei(p\x\)= I e ^-^ dt . (41) 



It is understood that all terms which vanishes with p — > are ignored. The leading non 
trivial contributions to F are 

Ei(/i \x\) = - 7 - Hp \x\) +0(p\x\), (42) 
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F 1 = e ^(-- + \x\{l-H^\x\)- 1 )\ ()[,,). 



8irr 2 \ fi 

F 2 = -^=(ln(/i|x|)+ 7 ) + 0(/i), 
87rV r 2 

Fg ^f(--\x\)(d-l) + 0(p), (43) 

2 2 2 2 2 

F = — (d - 2) + ^— + — ln^ W) - — b| Mil \x\) - — \x\(d-2 + 7), (44) 
8tt// ; 87rm 87rm VP 1 |; 8vr VP 1 |; 8tt 1 1 v /; ' v ; 

where 7 is Euler's constant and we set r 2 = m 2 .Here we used integrals for intermediate state 

for on-shell fermion 

J d 3 xexp(— ip ■ x) J d 3 rS(r 2 — m 2 ) exp(ir ■ x)f(r) = f(m), (45) 

d 3 xexp(— ip ■ x) I d 3 r exp(ir • x)5(r 2 — m 2 ) = — -. (46) 

J 4nx ' m 2 + p 2 

Exponentiation of F in(38) leads simple form 

exp(F) = exp(A - B \x\)(n \x\) D ~ c ^, (47) 

where 

A = —{d - 2) + 1?-, B = —(d-2 + -f),C = —,D = —. (48) 
8vr/i ; 8vrm' 8vr 8vr' 87rm V ; 



and we get the spectral function p[x) 

md 2 r 



mexp(- 


-m 


x\) 


47T 


x\ 





(2n) 2 ^/r 2 + m 



exp(ir-x), (49) 



mexp(- 


-m 


x\) 


47T 


X 





p{x) = — j— j exp(F(m,x)) 



mexp(-(mo + B) N) gxp(A)( | g | ) -CN +g> 

47T \x\ 



(50) 



2 

where m = \mo + ^{d — 2 + 7)1, B denotes the correction of mass,D is a coefficent of the 
bare Coulomb potential devided by m and C \x \ ln(/x |x|) can be understood as the position 
dependent self-energy as dynamical mass. Here exp(— m \x\)/A7r \x\ is a free scalar propagator 
with physical mass m.In Euclidean space we omitt the linear infrared divergent factor A. 
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B. Confining property 



Here we mention the confining property of the propagator Sf(x) in position space 



S F (x) = (- 



( 



m 
ij ■ d 



L)[ 



mexp(- 


-m 


x\) 


47T 


X 





+ l)p(x) 



D 



m 

e 2 e 2 
q _ _ 

8nm ' 87r 



(51) 
(52) 



Since is a dimensionless quantity and taking limit n — > with finite \x\ gives zero 

except for \x\ = oo.In section III we take the \x — > limit in momentum space. Thus here 



we fix /i = -umt of mass as m and see x dependence of the function (// |a;|) D ^'^^.The p(x) 
damps strongly at large x provided 



lim (fx\x\)- clxl = 0. 



(53) 



The profiles of the p(x) for various values of D > 1 are shown in Fig. 1. The effect of (/j, \x\) 

in position space is seen to decrease the value of the propagator at low energy and shown in 

Fig.2. 




Fig.l p(x) for 
m = \i = unit, D = 1, 1.5 



Fig.2 (p\x\)- Dm ^ for 
m = jj, = unit, D = 1 



In section III and IV we discuss dynamical mass, the renormalization constant, and bare 
mass in connection of each terms in F. 
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C. 0(e 2 ) propagator in Momentum space 



After angular integration we get the propagator 



= + (54) 

lit 

p{p) = -= \ rf|x|sin(v^2|a;|)exp(A- (m + B) \x\)(p\x\)- c W +D , (55) 

2iiJ-p 2 Jo 



where 



A = ^(d-2) + ^,B = f(d-2 + 1 ),C = f,D = ^—. (56) 
Slip Sirm 8n Sir Sttitl 

If we discuss the Euclidean or off-shell propagator we can omitt the linear infrared divergent 
part in A. In this case m denotes a physical mass 

e 2 

m= \m + — (d-2 + 7)|. (57) 
Here we show the propagator p(p) up to 0(e 2 ) 

<i 3 xexp(ip • x) / — — exp(ir • x)F(x) 

m 

= / d |x| sin( a/ — p 2 |x|) exp(— m \x\)[l + A — Cx\a{p \x\) + I? ln(/x |x|)] 

= [^4 + m(DI 1 -CI 2 )], (58) 
where I\ , J2 are the following integrals 



-7 ln((m 2 +p 2 )/p 2 ) ln((m - ^-p 2 )/(m + a/-P 2 )) 



m 2_|_p2 2(m 2 +p 2 ) m 2 +p 2 

sin (a/— p 2 |x|) exp(— m 



(59) 



f ^(v^k|)exp(-m|x|) |x| ^ w |x| 



-m 



: [ln(( m _ yf^)Hm + a/ 11 ^)) + ln((m 2 + p 2 )/^ 2 ) - 2(1 - 7)]. (60) 



(m 2 + p 2 ) 2 

From these expressions we see that the gauge dependent terms A,D,B, C,are wave func- 
tion renormalization constant and mass renormalization respectively.In this order the wave 
function renormalization constant 

„2 



Z {2) = i + A-- 
2 8vrm v ' 2 

is divergent at p 2 = —m 2 and p 2 = 00 



(7 + l ln((m 2 + p 2 )/p 2 ) + ln((m - yf^)/{m + a/ 3 ^)) (61) 
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D. Mass generation 



Since our model QED 3 is super renormalizable mass generation occurs perturbatively.In 
other words, mass changes as in the case of operator insert ion [8]. It is seen by the dimensional 
analysis of the propagator or self-energy S.For instance we take an example for the 0(e 2 ) 
self-energy 

£( 2 ) =1 . p A(p) + B(p), 

S (2) = e 2 J -0^1,S F (kh,D^ o (p - k), (62) 
7-pA(p)+B(p) 

Sf(p) = Aw + B*(py (63) 

whereE^) is given by setting A = 1,B = m for S^-By using contour integral we obtain 

-Lrr[( 7 -P)^) = A(p) = - P -^tzn\^)l (64) 

2jr 8np z mp m 

W ) ) = g(p)= (d "fi e2ro tan- 1 A. (65) 
2 Anp m 

Here we notice that in the limit B(p) p ^ = (d—2)e 2 / (47r) which is independent of m. From the 
above expressions we see that the high-energy behaviour of A(p) and B(p) are proportional 
to e 2 /p in this order with bare vertex. On the other hand,Dyson-Schwinger equation is non- 
linear 

d 3 k 1 B(k) 



B(p) = 2e 2 J 



(2tt) 3 (p- k) 2 k 2 + B(k) 
e 2 f°° kdkBik) ,p + k 



2tt 2 p J k 2 + Bikf^p-k^ (66) 
in the Landau gauge A = l.This yields B{p) — > m 3 /p 2 as p — > oo.We have the scalar part 
of the propagator 

1 e 4 
S E (p) = -^trS F {p) p ^ oc oc — , (67) 

which is the correct form of B(p) from dimensional analysis [8]. Here we return to our ap- 
proximation. The propagator in momentum space is 

Sf(p) = _ f ^■m(f) + ^ 2 ) = (TP + l)p(ph {els) 

J p 2 + ur + ie m 

P (p) = Jd* X eM-iVx)Wl (69) 
^ ^exp(- ffl |x|) (/jN) _ C | x | +fl {7Q) 
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From the above equations the conditions for mass generation is easily seen that 



< p(x) < Max, for 1 < D, (71) 
1 f°° 2 an(,/=j? \x\) 



pip = oo) = lim —f= / x 2 "*" vv f — p(x) = lim - — ir(x 2 p(x)) x=0 = 0, (72) 

P^OO y/- p 2 J Q \X\ P^OO yZ-pl 

p(p = 0) = J d 3 xpjx) = finite(~ 0.53 forD = l,m = n = l), (73) 

which are a prioli assumed to be satisfied in the numerical analysis of Dyson-Schwinger 
equation. 



III. ANALYSIS IN MOMENTUM SPACE 

To search the infrared behaviour we expand the propagator in the coupling constant 



e 2 and obtained the Fourier transform of p(x) [2]. In that case it is not enough to see the 
structure of infrared behaviour which can be compared to the well-known four dimensional 
QED. Instead we make Laplace transformation of (p \x\)~ c ^ : which leads the general spec- 
tral representation of the propagator in momentum space. After that we show the roles of 
Coulomb energy and position dependent mass. The former determines the dimension of the 
propagator and the latter acts to change mass. Let us begin to study the effect of position 
dependent mass (Self-energy), Coulomb energy in momentum space 

e 2 

exp(— \x\ M(x)) = exp(— — \x\ hn(p \x\), (74) 

57T 

eM Coulomb energy = * 

Similar discussion was given to study the effects of self-energy and bare potential in the 
stability of massless e + e~ composite in the lattice simulation [11]. The position space free 
propagator 

c i \ c a , ^exp(-m |x|) 

S F (x,m ) = -(«7 • d + m ) —— (76) 

47T \X\ 

is modified by these two terms which are related to dynamical mass and wave function 
renormalization.To see this let us think about position space propagator 

mexp(-mld). ,._*? M , , ,,^L_ 
Att \x\ 

It is easy to see that the probability of particles which are separated with each other in 
the long distance is supressed by the factor (p \x\)~ c ^,and the Coulomb energy modifies 
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the short distance behaviour from the bare 1/ \x\ to 1/ |x| 1_D .The effect of Coulomb en- 
ergy for the infrared behaviour of the free particle with mass m can be seen by its fourier 
transform[2,10] 



/ 



exp(— m 


M) 


Air 


x\ 





dxexp(-ip-x) — j— ; (p \x\) = p 



D D F(D + 1) sm((D + 1) arctan( v /^/m)) 



V^^ + m 2 )^)/ 2 
^{V^ - m)' 1 - near p 2 = -m 2 . (78) 



Above formula shows the structure in momentum space is modified for both infrared and 
ultraviolet regions. Usually constant D represents the coefficent of the leading infrared di- 
vergence for fixed mass in four dimension. Therefore Coulomb energy has the same effects in 
three dimension as in four dimension but change the ultraviolet behaviour since the coupling 
constant e 2 is not renormalized.Now we consider the role of M(x) as the dynamical mass at 



low momentum. First we define Fourier transform of the p(x); 

pip) = J exp(— ip ■ x)p(x)d 3 x. (79) 

Momentum dependence of p{p) for various values of D is shown numerically in Fig.3.We 
notice that 

D = - p{p) = -^- 2 , (80) 

D = 1 — > pip) = — — x log correction, (81) 

(p z + m z y 

which we can see from (70),(58)-(60).If we include (p \x\)~ c ^ term it is easy to see that 
the value of the proagator p{p = 0) decreases which is shown numerically in Fig. 3. In the 
previous paper we expand the propagator p{p) in powers of e 2 and see the logarithmic 
infrared divergence at p 2 = m 2 .In that case it is not clear the effect of position dependent 
mass. If we use Laplace transformation it is easily seen that (p \x\)~ c ^ acts as mass changing 
operator m — * m — s 

poo 

F(s)= / exp(-s\x\)(m\x\y c ^d\x\ , 
Jo 

POO 

(exp(-m* \x\)(m\x\)~ clxl = / exp(-(m* - s) \x\)F(s)ds, (82) 

Jo 



where we introduced 



m * =m + Chi(—), (83) 
m 
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to separate the cut-off p dependence from F(s) in the following way 

exp(— m \x)\ (/i Irrp ^' = exp(— m* \x\)(m \x\)~ c ^. 



(84) 



The Laplace transfom F(s) is shown in Fig. 4. 




Fig.3 p{p) for m — p — unit, D = Fig. 4 Laplace transform F(s) for 
0(upper), l(middle), 1.5(upper),p = (m \x\)~ Dm ^, 

n/20 m = 1,D = l,s = n/20 



We have the complete the expression of the propagator 

— - m/i D exp(-m* \x\) f°° hp , u 

p(x) = l— d / ex P( s \x\)F(s)ds, 



47T \X 



(85) 



D i T(D + 1) sm((D + 1) arctan( y^p^/(m* - s)))F(s)ds 

P\P) ~ m f J ' 



Pip = 0) 



V /= ^V / (P 2 + (s-m*) 2 ) 1+D 
°°2(s-m*)F(s)rfs 1 ^ +io ° 



(s-m*)< 



(86) 



2F(s)ds 

- — ~ 0.53 for m = = 1, D = 1. 

(s — m*) J 



(87) 



Here we study the cut off p — > limit. We see that m* diverges as ln(/x/m) in the limit 
p — > 0. In this case we can neglect the s dependence of the denominator of p(p) and using 



< 



s — m 



< 



s — m 



((p 2 + ( s ~ rn*) 2 ) 1+2 ° ip 2 + m* 2 ) +, 2 D 



poo poo 

/ F(s)ds — 1, / sF(s)ds = finite, 
Jo Jo 



(88) 
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we get the infrared behaviour 

- D) i 

\J—p 2 {p 2 + (s — m*) 2 )^~ m 



m 1+D e D T(l + D)sm((D + 1) arctan( ^/^/(s - m*)) // 

~ / — o, o — ; ^ t+d > 0? e ~~ ~ ~~ * °' ( 89 ) 



and for D = 1 

- ™Ml + !*(*)) - 0. (90) 
P{P> (p 2 + m 2 (l + ln(e)) 2 ) 2 1 J 

In Minkowski space we have by replacement p 2 — > — with infrared factor exp(A), 

2 2 

A = —(d-2) + —. (91) 
87r/i V 7 87rm V ; 

Since A is highly gauge dependent we do not have a definite answer in the limit /i — > 
except for the Yennie gauge d = 2. The infrared behaviour near p 2 = m 2 becomes 

_ T(D + l)e D sin((D + 1) arctan(p/m(l + D ln(e))) 
W ~ (-p 2 + m 2 {l + D\n{e)) 2 y+ D / 2 ' ' 



and for D — 1 



this leads to 



2m 3 e(l + ln(e)) , . 



p(p)« tt^tti - V ! pV™ !_ for£> = l,d = 2, (93) 
^ W m(ln(e)) 4 m (l n ( y/l - p 2 /m 2 )) 4 V ^ 

near p 2 = m 2 .In Euclidean space it is natural that in the long-distance the propagator 
vanishes in the limit \x — > O.In Minkowski space we find a vanishiment of the propagator 
as a/1 — p 2 /m 2 (ln(l — p 2 /m 2 ))~ 4 near p 2 = m 2 for the Yennie gauge at D — l,which is 
a consequence of confinement picture in our model. In comparison between our model and 
Dyson-Schwinger equation, D = 1 is prefered for dynamical mass generation.lt has been 
discussed that the propagator will have a branch point on the real p 2 axis associated with 
dynamical mass[12, 13]. Therefore we can say that the radiatively corrected fermion has not 
a simple structure as in QED 4 ,in QED 3 it has a superposition of different mass with cut and 
linear divergence associated with massless photon in the ordinary phase or it is expressed 
as superposition of dipoles with linear divergence in the condensed phase for finite cut-off 
/x, which will be shown in the next section. There is a possibility to remove infrared cut-off 
by including photon mass as Chern-Simon term or vacuum polarization of photon [3, 5, 7]. 
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IV. BARE MASS AND VACUUM EXPECTATION VALUE <^> 



In this section we examine the renormalization constant and bare mass and study the 
condition of vanishing bare mass based on spectral representation. The equation for the 
renormalization constant in terms of the spectral functions read 



Z 2 x (7 -p + mo) _ f 7 • ppi(uj 2 ) + p 2 (uJ 2 )du 
P ^oo p 2 — uiq + ie p^od 



t t i.c j p A — or + ie 

Instead we determine them directly by taking the high energy limit of Sp(p) 



(94) 



m Z 2 1 = m [ p 2 (uu 2 )duu 2 = - lim tr[p 2 S f{p)} (95) 

J 4 p^oo 

Z 2 l = [ Pl {u) 2 )du 2 = \ lim tr[ 7 • pS F (p)]. (96) 

J 4 p^oo 

In (86) the function sin((D + 1) arctan(A/— p 2 / (m* — s)) is expanded for large (— p 2 ), 



/V /= P^n AD + 1) N ,(-0 + 1) ,(D + l)(m*-s) 1 
sin((D + 1) arctan( V - ) = sin( V - ' tt ) - cos( V - ' tt)- £L= - + 0( 



mr — s 



From this we obtain for D < 1 



m 



Z' 1 = mp D T(D + 1) sin( (jP ^ 1)7r ) lim D 



Z 2 - 1 = ^r(D + l)sin( (jP + 1)7r ) lim 

Z v — +oo 



(0 < £>) 
m (0 = D) 

(0 < D) 

1 (0 = D) 



—p z 



(97) 

(98) 
(99) 



For D — 1 case first term in the 1 / \/—p 2 expansion vanishes and 

m^Z 2 x = lim ^ / (m* - s)F(s)ds = 0, 
p^oo — y 

., . (m* - s)F(s)ds = 0. 



Z 2 1 = lim 



2/j 



(100) 
(101) 



p^oo — p^ 

This means that propagator in the high energy limit has no part which is proportional 
to the free one and it shows confinement. Usually mass is a parameter which appears in the 
Lagrangean.For example chiral symmetry is defined for the bare quantity.In ref[9] the relation 
between bare mass and renormalized mass of the fermion propagator in QED is discussed 
based on renormalization group equation with the assumption of ultraviolet stagnant point 
and shown that the bare mass vanishes in the high energy limit even if we start from the 
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finite bare mass in the theory.lt suggests that symmetry properties can be discussed in terms 
of renormalized quantities. In QCD bare mass vanishes in the short distance by asymptotic 
freedom. And the dynamical mass vanishes too [8]. Since our model is super-reomalizable bare 
mass m does not vanish. In our apprximation this problem is understood that at short 
distance propagator in position space tends to 



p( x ) = 7, , U| u exp(-B\x\+D\n(fi\x\)-C\x\\n(fi\x\)) x ^ (102) 

47T \x\ 

-> |x| D-1 (Ai|a;|)- c|ie| , (103) 



where we have p(0) = finite at D = 1 case which is independent of the bare mass mo-Thus 
we have a same effect as vanishing bare mass in four dimensional model. Of course we have 
a dynamical mass generation which is m = |f^(d — 2 + 7)| + M(x) for m = in our 
approximation. There is a chiral symmetry at short distance where the bare or dynamical 
mass vanishes in momentum space but its breaking must be discussed in terms of the values 
of the order parameter. Therefore it is interesting to study the possibility of pair condensation 
in our approximation. The vacuum expectation value of pair condensate is evaluated 



exp(- 


-m 


x\) 




X 





») = -trS F (x) = -2m/i D lim ^ K , ^^ (// \x\^ c ^ +D 

^0 (D>1)^ 

finite (D — 1) , 
v oo (D<1) J 



(104) 



for finite cut-off //.For D < 1 case the p{x) is divergent at x = and it looks like a 



spike, for D = l,p(0) = finite and it looks like a wave packet at finite range where long 



range correlation appears, and for D > l,p(0) = but it does not damp fast for large x.ln 
the weak coupling limit we obtain Z 2 = l,m = m and (ipip) = oo.If we introduce chiral 
symmetry as a global C/(2n),it breaks dynamically into SU(n) x SU(n) x U(l) x U(l) as in 
QCD[8,11] for D — 1 for finite infrared cut-off.Our model may be applicable to relativistic 
model of super fluidity in three dimension. Usually we do not find the critical coupling D = 1 
in the analysis of the Dyson- Schwinger equation in momentum space where only continuum 
contributions are taken into account and we do not define physical mass. 
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V. SUMMARY 



Infrared behaviour of the propagator has been examined in the Bolch-Nordsieck-like ap- 
proximation to QED 3 . Using LSZ reduction formula for the one photon matrix element we 
find two kinds of logarithmic infrared divergent terms. These are position dependent mass 
and Coulomb energy both of which are gauge independent in our approximation. Confining 
property is shown by position dependent mass at large distance where propagator vanishes 
faster than the sloution of Dyson-Schwinger equation converted into position space with 
some approximation. Renormalization constant Z^ 1 vanishes for arbitrary coupling D in our 
approximation, which implies confinement. Momentum space propagator is expressed as dis- 
persion integral and it shows us that the effects of dynamical mass can be written as a 
superposition of the propagator with different masses. In our approximation there seems to 
be a critical coupling constant for vacuum expectation value (ipip) which is independent of 
the bare mass. Our model shows that the confining property of charged fermion and dynam- 
ical mass generation are realized in QED 3 (at critical coupling constant D = 1) as in QCD 
in four dimension with finite infrared cut-off which may be supplied by some dynamical 
origin of the photon mass beyond the quenched approximation . 
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VII. A. GAUGE DEPENDENCE 

In ref[2,5,9] gauge dependece of the propagator has been discussed. In order to understand 
the gauge dependence of the propagator away from the threshold, we will once examine the 
gauge covariance relation. Writing D /xl/ (k) = D ( $(k) + k il k v M{k) ,we see that the fermion 
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proagator varies with the gauge according to [4], 

S F (x, K) = exp(-iKy/a?)S F (x, 0), (105) 

—— — (exp(-ik -x)-l) = -iKVx*, (106) 
[2n) 6 k 4 

where K = e 2 <i/87r.Thus the gauge transformation property is violated by linear infrared 
divergences near the threshold. Here we adopt the subtracion of the point k = in the 
integral. Free propagator with mass u> equals to 

a(.°)f \ c flj xexp(-^ v /Z a^) nrV7 , 
Sy(x) = -(i^-d + u) — 4?r y-^ — ' ( 107 ) 

and its Fourier transform is 

SP(p) = j d 3 xexp( ip ■ x)SP(x) = ^ p2 P _ + J - (108) 

In our approximation only mass terms changes under the gauge transformation and re- 
maining terms in F are gauge invariant. In general gauge invariant part is approximation 
dependent. 

Finally we show the general spin dependent spectral function in the 0(e 2 ) for d gauge by 
evaluating 

F v = -£-2 tr (-Y-rF),F 9 = ±tr(F), 

and we have 

F = (^Z+ l)[f(-(d - 2 + 7 ) \x\ - f \x\ Info |x|) + /-(Info \x\) + 7 )] 
m 87r 87r sum 

i_LptJl. (109) 

m Sirm* 2 \x\ 

The first term is the same as the previous one. The second term is not significant in the 
infrared. This term is also obtained by gauge transformation from the Landau gauge where 
Pi = p 2 is assumed [7]. Since 

7 ■ r + m 
2m 

is a projection operator it is easy to show the following. In Euclidean space if we exponentiate 
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we have 



S F (x) = - I !^— 2 exp(.r ■ x) (1 + ^) exp((^ + 1)F) 
7 v r 2 + m 2 mm 

f mrdr f • ^ . , r '7 • w-^\\ 

= ~ / / 2 , 2 exp(zr -g)(l + )exp(-)(cosh(-) + sinh(— )) 

7 Vr 2 + m 2 ml 2 m 2 

mrdr r • 7 . . . . 

exp(ir-x)(l + -)exp(F), (110) 



Vr 2 + m 2 m 
S F (p) = J d 3 xexp(-ip-x)S F ( y x). (Ill) 
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